Minimization of Ohmic losses for domain wall motion in a ferromagnetic nanowire 
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We study current-induced domain-wall motion in a narrow ferromagnetic wire. We propose a way 
to move domain walls with a resonant time-dependent current which dramatically decreases the 
Ohmic losses in the wire and allows to drive the domain wall with higher speed without burning the 
wire. For any domain wall velocity we find the time-dependence of the current needed to minimize 
the Ohmic losses. Below a critical domain- wall velocity specified by the parameters of the wire the 
minimal Ohmic losses are achieved by dc current. Furthermore, we identify the wire parameters for 
which the losses reduction from its dc value is the most dramatic. 
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Introduction. In recent years there has been intense 
interest in applications of domain wall (DW) motion in 
ferromagnetic nanowires [ll, This interest is mostly 
based on the possibility to store and exchange informa- 
tion by means of moving domain walls which separate 
the regions of magnetization parallel and anti-parallel to 
the wire. These regions with parallel and anti-parallel 
magnetization can be thought of as two bits, zero and 
one, of binary information storage. 

DWs can be moved by a magnetic field @, Q or electric 
current [l|, 01 ■ For technological applications the current 
driving is preferred as magnetic field is difficult to apply 
locally to small wires. Thus, in this Letter we consider 
the current-driven DW devices. To achieve their high- 
est performance it is important to minimize the losses 
on Joule heating in the wire, which are due to the resis- 
tance of the wire itself and the entire circuit. They are 
proportional to the time-averaged current square, (J'^). 
Their minimization has a twofold advantage. First, one 
can increase the maximum current which still does not 
destroy the wire by excessive heating and therefore move 
the DWs with a higher velocity, since the DW velocity 
increases with the applied current. Second, it creates the 
most energy efficient memory devices and also increases 
their reliability. 

To achieve these goals we propose to utilize a "res- 
onant" time-dependent current, which allows to gain a 
significant reduction of Ohmic losses. We show that all 
thin wires can be characterized by three parameters ob- 
tained from dc-driven DW motion experiments: critical 
current jc, drift velocity at the critical current Vc, and 
material dependent parameter a > 0, which in particular 
depends on Gilbert damping a and non-adiabatic spin 
torque constant /3. The parameter a is just a ratio of the 
slopes of the drift velocity Vd{J) at large and small dc- 
currents, see the upper inset of Fig. [TJ Our main results 
are summarized in Fig. [5] We find the minimal power 
(J^) needed to drive a DW with drift velocity Vd- Fig- 
ures [2] (a) and (b) show the dependence of power ( J^) 
on Vd for the optimal time-dependent current - red solid 
curves, and for dc current - black dashed curves, for two 
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FIG. 1: (color online) A sketch of a current driven domain 
wall in the ferromagnetic wire. The upper inset shows the 
dependence of drift velocity Vd of DW on dc current J for 
B > and B < 0, see Eq. (f2b)) . The slope at J < jc is given 
by A, while at J 3> jc it is A + S. The lower inset shows 
the power of Ohmic losses Pdc{Vd/Vc) = /jc for dc current. 
For B < the power has a discontinuity at Vd/Vc ~ 1. 



cases: (a) a < 1 and (b) a > 1. In Fig. [5] (a) the minimal 
power is given by dc current for Vd < Vc, but above Vc 
there is a significant reduction in the heating power com- 
pared to dc current. Fig. [21(b) shows that the power (J^) 
is reduced in comparison with the dc case for Vd > VcV^c- 
The (dimensionless) resonant critical velocity v^c < 1 and 
can be extracted from the dc-current measurements. For 
Permalloy using Q we estimated a ~ 0.5, see Fig. [5] (a), 
where for Vd ^ Vc the power is less than 50% of that for 
the dc current. 

Figs, m (c) and (d) show the limiting cases of a <C 1 
(c) and a ^ 1 (d). We note that for small a and /3, 
a sa a/13. If a < 1 (/3 > a), Fig. [H (c), for dc current 
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FIG. 2: (color online) Minimal power of Ohmic losses p — 
{J^)/jc as a function of drift velocity Vd shown by solid line 
for (a) a = 0.5 (b) a — 2. The dashed line depicts p for dc 
current. A sketch of {J^){Vd) shown by solid line in (c) for 
/3 » a (a < 1) and (d) for /3 < a (a > 1). 



the excessive heating power ^ essentially limits the 
highest achievable drift velocity Vd by Vc, whereas the 
resonant ac-current can move DWs with much higher Vd 
(and still rather low power). In the opposite case a 3> 
1, (/3 <C a), Fig. [5] (d), the power saving starts to be 
considerable at very small velocity Vd- If /3 = the dc- 
current power is finite even at Vd — 0, while for the 
resonant ac-current the power linearly approaches zero 
at small Vd- Therefore, our approach gives a dramatic 
power reduction even in the least favorable cases /3 ^ a 
and /? 3> a, thus opening new doors for using materials 
with much wider range of /? for fast DW motion. 

Model. DW in a ferromagnetic wire can be modeled 
by a Hamiltonian which contains exchange and dipolar 
interactions. In a thin wire, the latter can be approxi- 
mated by two anisotropics: along the wire (A) and trans- 
verse to it (K). A sketch of a wire with a DW of width 
A is shown in Fig. [T] The dynamics of magnetization S 
in a wire is described by Landau-Lifshitz-Gilbert (LLG) 
equation with the current J 0, 0) ; 

S = S X He - jaS + /3JS X aS + aS X S, (1) 

where Hg = 6H/dS is the effective magnetic field given 
by the Hamiltonian T-L of the system, a is Gilbert damp- 
ing constant, l3 is non-adiabatic spin torque constant, and 
d = d/dz. Furthermore, it can be shown [Sf that in a thin 
wire the DW is a rigid spin texture for not too strong ap- 
plied currents and its dynamics can be described in terms 
of only two collective coordinates (corresponding to the 
two softest modes of the DW motion), namely, the posi- 
tion of the DW along the wire zq and the rotation angle 
4> of the magnetization in the DW around the wire axis. 

To describe the DW dynamics we need to find the equa- 
tions of motion. For the two softest modes of the DW, 
zo(0 s-^d (j){t), they can be found as an expansion in small 



current J up to a linear in J order. Due to the transla- 
tional invariance zq and 4> cannot depend on zq. In addi- 
tion, to the first order in small transverse anisotropy K, 
<j) and io are proportional to the first harmonic sin(20). 
Then the most general equations of DW motion are 



C[J-jeSin(2(/.)], (2a) 
= AJ-KB[J- jcsin(20)], (2b) 



where J{t) is, in general, a time-dependent current whose 
frequency is not too high to create spin waves and other 
excitations in the wire. Coefficients A, B, C, and critical 
current jc can be calculated for a particular model 14 1 
in terms of a, /3 and other microscopic parameters by 
means of deriving Eqs. ([2]) from LLG equation ([1]). How- 
ever, we emphasize that Eqs. ([2]), with coefficients A, 
B, C , and jc determined directly from dc-current experi- 
ment for each particular wire, have more general validity 
than just being derived from LLG, e.g., due to the com- 
plicated influence of disorder and internal DW dynamics 
(9[ . Namely, the value of jc is deflned as the endpoint 
of the linear regime of the time-averaged (drift) velocity 
Vd = {zo{J)), see the upper inset of Fig. [TJ The lin- 
ear slope of Vd{J) below jc determines constant A. The 
slope of Vd{J) at large J gives A + B. Constant C one 
can obtain, e.g., from the measurements of the DW elec- 
tromotive force III for dc current. 

DC current. For the dc current applied to the wire 
the DW dynamics governed by Eqs. can be obtained 
explicitly lS\ . For J < jc and A ^ the DW moves along 
the wire but does not rotate around its axis. It only 
tilts on angle 00 fi'om the transverse- anisotropy easy axis 
{y axis) given by condition sin(2(/)o) = J/jc- The drift 
velocity is given by Vd = AJ, see Eq. ((2b|) . At J = jc the 
magnetization angle becomes perpendicular to the easy 
axis, 00 = Tr/2. For J > jc the DW both moves and 
rotates, and Eqs. ([2]) give Vd — AJ 



Jc 



The influence of the spin structure on the current 
is negligible. The largest losses in the system are the 
Ohmic losses of the current. The power of Ohmic losses 
is proportional to J^. Therefore, at J < jc the cur- 
rent is J = Vd/A and the power of Ohmic losses is 
'Pdc — = VljA^. It is instructive to introduce the 
dimensionless variables for time, drift velocity, current, 
and power. Using Vc — Ajc — K/S. we flnd [15] 

T^cjct, vd^Vd/Vc. j^J/jc, P = r/jl (3) 

Using Eq. ([3]) we flnd pdc = v'^ for < 1. 

For currents above jc the dimensionless power pdc is 
given in terms of dimensionless drift velocity Vd — j + 
(BlA)\/j'^ — 1 as Pdc(wd) = i^, see the lower inset of 
Fig. [1] Thus, it is quadratic in Vd, and at 3> 1 it 
approaches pdc = J^v\l{B -f ^ B/{B + A). For 
B > Q right above = 1 , it is approximated by pdc = 1 + 
{A/ B)'^{vd~iy' . For B <Q the power has a discontinuity 
a,t Vd = 1. 
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Minimization of Ohmic losses by time- dependent cur- 
rent. In this part we minimize the Ohmic losses while 
keeping the DW moving with a given drift (average) ve- 
locity. Equations of motion ^ are correct even when 
the current depends on time. In general, the DW motion 
has some period T and current j(r) must be a periodic 
function with the same T to minimize the Ohmic losses. 

In the following it is more convenient to measure the 
angle from the hard axis instead of easy axis and to scale 
it by factor of 2, so that 2(j) — O — n/2. Also, we introduce 
the ratio of slopes of Vd{J) at large and small currents 
a — {A + B)/A. Then using Eq. (l2a| the dimensionless 
current becomes 



3{T) = e/2-cose, 



(4) 



where 9 = 89/ dr. Averaging Eq. (j2bl) over dimensionless 
period T we find 



Vd 



(0) 



(COSf 



(5) 



where (...) = ^ Jq^ . . . dr is the time averaging. 

To minimize the power of Ohmic losses p averaged over 
time we need to find the minimum of {P{t)) at fixed Vd 
given by Eq. ([5]), 




-2p[^9^cos9-Vd'^ 



(6) 



Here to account for the constraint given by Eq. ([S]) we 
used a Lagrange multiplier 2p, with p being an arbi- 
trary dimensionless constant. Note that the cross term 
~ / ^ cos 9 dr' and the term ^ J 9 dr' can be dropped for 
the minimization procedure as they are full derivatives. 

Power ^ can be considered as an effective action for a 
hypothetical particle of mass 1/2 in a periodic potential 
field, and its minimization leads to the equation of motion 



dU 



cos2 6'-2pcos6'. (7) 



It can be reduced to the first order differential equation 

9 = ±2^d~Ui9,p), (8) 

where d is an arbitrary integration constant. Note that 
changing p — > —pinU of Eq. (jT)) is equivalent to changing 
9 ^ n -\- 9, so below we consider only positive p. The 
potential has a minimum at 6* = with Umm = 2/9—1 
for any p > 0. For p < 1 it has also minimum at 9 — ±7r 
with [/(±7r) = — 2p— 1 and the maximum at cos 9p ~ — p 
with U {i9p) = p^. For p > 1 it has maximum a.t 9 — zLtt 
with C/(±7r) = 2p- 1. 

According to Eq. (jS]) there are two different regimes: 
i) the rocking regime where d < max[U(0, p)] in which 
case 9 is bounded, and the particle oscillates in potential 
well U{9), see Fig. [3] and ii) the rotational regime where 




FIG. 3: (color online) Potential U{6) in which a "particle" is 
moving in the rocking (pendulum-like) and rotational regimes. 



d > max[U(^, p)] in which case the magnetization in the 
DW rotates. Below we consider these regimes separately. 

Rocking regime. In this regime the motion of 9 mimics 
pendulum motion. The particle rocks between the two 
turning points —9q and given by the condition d = 
U{±9q,p). At these points 9 = 0. Since is a bounded 
function (9) = and the averaged velocity becomes Vd = 
— {cos 9). The averaging is done over a period of one 
complete oscillation, 



T 



d9 



-00 y/d^U{9,p) 



(9) 



and according to Eq. ^ the power is given by 

p= {9^)/4+ {cos^9). (10) 

Most generally 9 depends on time. For any 9{t), how- 
ever, (cos^0) > (cos0)^. Then from Eq. (fTU|) follows 
p> (^2)/4-i- (cos6l)2 = {9^)/4 + vj > vl=pdc. Thus, in 
the bounded regime the power of Ohmic losses is minimal 
for dc current and is given hy p = v^. 

Rotational regime. Next we study the case when 
d > max[U(6', p)], so that angle 9 is unbounded. It corre- 
sponds to the rotational motion of the transverse to the 
wire component of the DW magnetization. Note that in 
the rotational regime the term in Eq. ([5]) with (9) should 
be kept because 9 is not bounded. The time it takes for 
9 to make a full rotation from — tt to tt defines the period 
T. Then the period, drift velocity, and power, according 
to Eq. pO|) . are given by 



1 



T = - 



d9 



I'd = 777 



2.U y/d-U{9,py 

1 r COS 9 d9 
'2T J_^^d-U{9,py 
1 r d- U{9, p) + cos^ 

2T .Lr^/T=www 



na 
T 



d9. 



(11a) 
(lib) 
(11c) 



This system of equations, after minimizing the power p 
with respect to both d and p at fixed Vd, gives p{vd). One 
can either directly perform a numerical minimization of 
Eq. (|lip or alternatively try to find the minimization con- 
dition for p analytically. We have followed both routes. 
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FIG. 4; (color online) Current j as a function of time r at 
velocities = 0.5 (dashed line), «d = 1.5 (dot-dashed line), 
and Vd = 4.5 (solid line) in the rotational regime for a = 2. 



The minimization of Eos. (fTT]) infers that dp/dp\v^ — 
from which one can find \m 



y^d-Ui0,p)d9 = 27rap. 



(12) 



This equation gives the relationship between d and p. 
Solving it together with Eq. (|llbp one finds d and p in 
terms of Vd- They are then substituted into p = 2pvd — d 
which follows from Eqs. (|llc|) and ([T2|) . The motion is 
unbounded when d > max[U(0, p)], see Eq. ([5]), which 
leads to d > p^ for p < 1 and d > 2/9 — 1 for p > 1. 

The results for the minimal power of Ohmic losses 
p{vd) are presented in Fig. [21 For a < 1, see e.g. Fig. [2] 
(a), at Vd < Wrc = 1 the minimal power p coincides with 
the one given by dc current, whereas at > 1 it is signif- 
icantly lower than pdc • Immediately above = 1 we find 
that there is a range of Vd where p — 1 + 2po{vd — 1) with 
po{a) > 1 given by Eq. ([12]) with d = 2p - 1. Therefore, 
p is linear in Vd right above Wrc = 1. 



For a > 1, see e.g. Fig. [2] (b), we show [l2| that there 
is a critical velocity Vrc < 1, such that at Vd < Wrc the 



power of Ohmic losses is p 



Pdc- Above Vrc one can 



minimize the Ohmic losses by moving DW with resonant 
current pulses. Right above Vi-c there is a certain range 
of Vd where d p^ , and therefore we find p = 2pQVd — Pq 



with po{a) < 1 given by Eq. (|T2|) with d — p^ . The crit- 
ical velocity is found as v-^c — Po(a)- For a ^ 1 (corre- 
sponding to non-adiabatic spin transfer torque coefficient 
/3 ^ a, c.f. Eq. ([T])) we find Wrc — 2/{na) and therefore 
for Vd ^ Vrc we obtain p = Avd/{Tra). 

We show that at large Vd the minimal power is always 
smaller than p^^- Note that for d 1 Eq. (IT^ gives 
d — a^p^. Using it we find that the difference between 
them approaches pdc — P = (1 — 1/^)^/2 at ^ 1. 

Optimal current. For Vd < Vj-c the optimal current 
coincides with the dc current. Above v^c the resonant 
current j{t) is plotted in Fig. H] for different velocities Vd 
in the case a = 2. At small Vd the current is given by 
j(r) = — 2cos(0(t)) — Vrc for cos(6'(r)) < —v^-c and by 



j(T) = Vie for cos(6'(r)) > — Wic- At Vd ^ 1 the current 
is approximated by j Vd/a + [(1 — a) /a] cosO. 

In general, at Vd > Wrc the current's maximum jmax in- 
creases from 2 — Vi-c at small enough ^ 1 up to jmax ~ 
Wd/a at 3> 1. The current's minimum increases mono- 
tonically from small positive values jmin = Vrc at 1 
up to jmin = jmax - 2|1 - a|/a at Vd > 1 . At < 1 (for 
a > 1) the time between the current picks decreases with 
increasing velocity as T ~ (tta — 2 arcsin'i;i.c)/(W(i — Wrc), 
whereas the pick's width is given by ~ 1.3/ (1 — Wrc), 
which is independent of Vd. Therefore, at small Vd — frc 
the picks are widely separated, then as Vd increases the 
time between the picks decreases. At S> 1 the op- 
timal current has a large constant component, which is 
close to but smaller than the dc current for the same Vd, 
and has small-amplitude ac modulations with a period 
T « na/vd on top of it. 

Summary. We have studied the current driven DW dy- 
namics in thin ferromagnetic wires. We have found the 
ultimate lower bound for the Ohmic losses in the wire for 
any DW drift velocity Vd. The explicit time-dependence 
of current, see Fig. [51 has been found which minimizes 
the Ohmic losses. We have shown that the use of these 
specific current pulses instead of applying dc current can 
help to significantly reduce heating of the wire for any 
Vd. Even in the limiting cases of the systems with weak 
(/3 <C a) or strong {(3 3> a) non-adiabatic spin transfer 
torque, where the power of Ohmic losses is high for dc 
currents, the optimized ac current gives significant reduc- 
tion in heating power thus greatly expanding the range of 
materials which can be used for spintronic devices [l|, 0] ■ 

We are grateful to J. Sinova for valuable discussions. 
This work was supported by the NSF Grant No. 0757992 
and Welch Foimdation (A-1678). 
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[15] It can be shown that C ^ B ^ a — P 0]. In the special 
case of a = P, we find C — B — and one cannot use 
dimensionless variables (|3]). The DW dynamics in this 
case is trivial [13] . The DW does not rotate = and 
moves with the velocity given by current io = J- 



Supplementary material for "Minimization of Ohmic 
losses for domain wall motion in a ferromagnetic 
nanowire" 

MINIMIZATION PROCEDURE 



As described in the main part of the Letter we study 
a domain-wall dynamics under the influence of a time- 
dependent current. The equations of motion for the do- 
main wall in the thin ferromagnetic wire take the form 
of Eqs. (2) of the main part of this Letter. To obtain 
the results of the main part, we use the dimensionless 
variables introduced in Eq. (3) of the main part. 

In order to minimize the power of Ohmic losses we need 
to find a minimum of the average j'^, 

P={fir)) (13) 

at a fixed drift velocity v^. From Eqs. (2) of the main 
part of the Letter we find 



Jir) = f- 




cos 9, 
cos 9, 
(cos 9), 




(14) 
(15) 
(16) 

(17) 



where 9 = 89 /dr and the averaging is performed over 
the dimensionless period T of the magnetization oscilla- 
tions. To find the minimum of power at fixed drift veloc- 
ity we introduce a Lagrange multiplier 2p and minimize 
the functional 

-cos 9] \-2p(^^{9)-{cos9)-Vd), (18) 



where (...) = T ^ ... dr is the averaging over time. 
Then it follows 

p=^^ (^^9^ +cos^ 9 + 2p cos 9^ dr, (19) 

where we dropped all the terms which are full deriva- 
tives and therefore do not give any contribution to the 
minimized power. 

The minimization of Eq. (1191) gives the following equa- 
tion of motion 

^9=-deUi9,p), (20) 
U{9,p) = -cos"^ 9 -2pcos9. (21) 
Its solution is given by 



9 = ±2^/d-U{9,p) (22) 
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where d is an arbitrary constant of integration. We note, 
that changing p — > — p is equivalent to changing — 
Tr+9, so we only need to consider positive p. The extrema 
of potential U are 



0, 

±7r, 



cos 9o 



U{0,p) = 
U{±7r,p) . 

u{±ep,p) 



l-2p, 



(23) 



where the last extremuni exists only for p < 1. 

We see that there are two different cases: i) rocking - 
when d is smaller then the maximum of U{9,p), and ii) 
rotating - when d is larger then the maximum of U{9, p). 
We consider them separately. 



ROCKING REGIME 

In this case there is an angle 9q such that d= U{9o,p). 
Then we have the oscillating motion between —9q and 
and back. The period of one complete oscillation, as well 
as the drift velocity and power are given by 



T = 



Vd 



d9 



-oo ^U{9o,p)-U{9,p)' 



- cosf = -- 



cos 9 d9 



T 7_e„ ^U{9„,p)-U{9,p) 



-9'- + cos^ 



l_ Z""" U{9Q,p)-U{9,p) + cos^i 
T J.eo ^U{9o,p)-U{9,p) 



■d9. 



(24) 



,(25) 



(26) 



dc current 

In the case of dc current 9 = const equation (|22|) gives 
9 = 9o, then ([20|) requires that is an extremum of 
U{9,p), there are three possibilities 9 — {0,TT,9p}, for 
them we get Vd = {—1, 1, p} = j, and p = {1, 1, p^}. So 
we see, that p — v'^. 



Thus, we see that in the rocking regime the dc current 
minimizes the Ohmic losses. 



ROTATING REGIME 

In this section we study the rotational regime. In this 
case 9 makes a full rotation from — tt to tt. The period of 
one complete rotation, as well as the drift velocity and 
power are given by Eqs. (11) in the main part of this 
Letter, i.e. 



2 J-^ ^d-U{9,py 



td9 



1 • 



^ ' cos^( 



T 2T ^d~U{9,p) 
1 r d-U{9,p)+cos^9 
2tJ_^ ^d-U{9,p) 



(30) 
, (31) 

d9, 
(32) 



where a = {A + B)/A. 



dc current 

In this case j = const and according to Eq. ([M]) we 
find 

1- , f'^'' d9 TT 

-9-cos9^j, T= ^ , 

2 70 J+COS0 

(cos 9) =-j + \ (0) , Vd = ^ (9) + J. (33) 

Since (9) ~ 2tt/T we obtain 
a- 1 



Vd 



Vf-^ + J, Pdc=f, j>l- (34) 



In particular for j 3> 1 we find 



ac current 



Now we take 9 changing with time. From the condition 



((008 6*) - cos9f'^ > 0, 
it follows that 

{coa^9) > (cos 6*) 2. 
Then the power p satisfies the condition 

p = ^-(^9^) + {cosH)>l{9^) + {cos9f 



''^+vj>vj 



(27) 
(28) 



(29) 



Pdc 



(1 + a) 



4 , a - 1 
vi ' 



a + 1 



ac current 



The same trick as in Sec. □ gives 



> 



1 + a^ 



Vl + a? ■ a 



(35) 



and we can only conclude that the power is no smaller 
than quadratic in Vd- 
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Derivation of minimization condition, Eq. (12) 

In this section we derive the minimization condition, 
Eq. (12) of the main part of the Letter. We rewrite 
Eqs. (11) of the main part in the foUowing way 



~ 2 ^d-U{9,p) 



TVd = TTfl - - 

p = 2pvd - d + 



1 r cos Ode 
1 ^ 



(36) 
(37) 



T 



y/d~U{9,p) - ap dO. (38) 



The minimization of p means that dp/dp\vd = 0, and 
we find 



dp 




dT 


1 r - 


dp 




" 'dp 


r2 / 

Vd "'-T 



^d-U{e,p)-ap 



dO^Q 
(39) 
I 



There are two possibihties to satisfy this condition 

= 0, or (40) 
dO = 0. (41) 



dT 
~d~p 



^yd-U{e,p)-ap 



Note that there is also a possibihty T — oo but it corre- 
sponds to a dc-current case. 

First we consider the possibihty given by Eq. (|40l) . Dif- 
ferentiating Eqs. (|36p and ((37|) with respect to p at fixed 
Vd and using Eq. (|40p . we obtain 



op 
op 



dO 



^ [d-C/(0,p)]3/2 

cos6'd6' 



^dO 



, [d-C/(0,p)]3/2 
COS^0d0 



Combining these two equations, we find 



cos d0 



n 2 



„ [d-C/(0,p)]3/2 



^6* 



cos2 61' dO' 



, [d-C/(0,p)]3/2 [d-C/(0',p)]3/2 
I 



(42) 



Note that this is a standard Bunyakovsky (Cauchy ~ 
Schwarz) inequahty and it never becomes an equality 
except for = Q. Therefore, we conclude that in the 
rotational regime the minimization condition is given by 
Eq. (liTI) . Then, the system of equations which we need 
to solve is 



y/d-Uie,p)-ap 



dO = 0, 



VdT = TTfl - - 

p = 2pvd - d. 



cos 61 ^6* 
. ^d-U{0,py 



(43) 

(44) 
(45) 



where the second equation is just Eq. p7l) with T given 
by Eq. ((36l) . Eq. (|43| provides the correspondence be- 
tween parameters p and d. Solving together the system of 
equations (l43|) and (j44|) yields /o(i;d) and d(t;d)- They are 
then substituted into Eq. (j45|) to find p{vd)- In general, 
the system of equations (^5]) and dH]) has to be solved 
numerically. However, in the limiting cases of small and 
large drift velocities Vd the result for p{vd) can be ob- 
tained analytically. Below we solve these limiting cases. 



Large Vd 



Using it we find that the difference Pdc — p = (a— l)^/2a^ 
&t Vd^ 1. Thus, at large Vd the minimal power is always 
smaller than pdc- 



Small Vd 

We recall that in the rotational regime the motion is 
unbounded and therefore d > max[U(6', p)], see Eq. (|22p . 
This leads to d > for p < 1 and d > 2p — 1 for p > 1. 
We show that it is possible to find analytical results for 
p{vd) when d = -t-e or d = 2p — 1 -f e and e ^ 1. Below 
we consider these two cases separately. 



< d < 1 



For the case d = p^ + € with e ^ 1, we find 



d0 



+ (p + cos0)2 



T(e, p)vd = TTfl - 



1 



cos d0 



2 J-^ (p-f cos6')2 



(46) 
(47) 



First we consider the case of large Vd which corresponds 
to large parameters d. At d >• 1 Eq. P5]) gives d = a^p^. 



In particular if p = 0, one can find that T = ln(16/e), 
Vd — 7ra/ ln(16/e), and p = 4/ ln(16/e) = Avd/ {ira). 
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From Eqs. (ge]) and (gT]) we find 

[Vd - p)T = ira - - / , , , , == (48) 



(49) 



Also, substitution of d = + e into Eq. (|43p yields 



/ ^e + {p + cosefde = 
Jq 



nap, 



So far equations and are exact. Note that the 
right-hand side of both equations (|48l) and (l49l) converge 
even if we set e = 0. Therefore, setting e = in them, 
we obtain 

[Vd - pjT = an - - I - — ——^dO, 



2 \p + cost 



and finally 



where cosOp 
rewritten as 



\p + cos 9\d9 = nap, 



(vd- p)T = n{a + l)-29p, 
2p0p — np + 2 sin 6p = nap, 



(50) 
(51) 



-p. This system of equations can be 



T = 
tan^ 



n{a + 1) - 26'p 



Vd 



(1 + a), 



(52) 
(53) 



and gives the period T and parameter p (since p = 
— cosOp). Since we look only for a solution of Eq. ([55]) in 
the range n/2 < 9p < n, we note that this solution exists 
only for a > 1. Also, as T must be positive (otherwise 
there is no solution for e) we see that Vd > p- 

Thus, we conclude that at a < 1 the minimal power for 
f < 1 is always achieved by dc current. At a > 1, based 



on Eqs. (1521) and (|53p . we find that p = v^-c for Vd > frc 
(as long as e is still small), and p = Vd for < Vd < Wrc- 
Therefore, for Vd < Urc, P = ^d and period T — go which 
corresponds to the dc current case. For Vd > Wrci 



p = 2vrcVd 



(54) 



and this power is minimized by the resonant ac current 
with the period between picks given by 



2 tan 9y^^ 

Vd - Wrc 



(55) 



In the limiting case of a ^ 1 we obtain p w 2 /an <^ 1, 
and for Vd^ p find p « Avd/an with T « an jvd- 



d > 1 

Now we consider the case ofc? = 2/3— 1 + e with e <C 1 
which corresponds to Vd immediately above v^c — 1 for 
a < 1. Neglecting very small e the power becomes 

p(vd) = \ + 2p^{vd-\). (56) 
Note that the power p is linear in Vd right above v-s-^ — X. 
Also for small a one can see that at ii^ > 1 the power is 
significantly lower than pdc- 

The parameter po in Eq. ()56|) can be found from 
Eq. gS]) with d = 2p-\. 



V(l + cos0)(2po + cos(? — 1) + e o?6' = Trapo- (57) 



In the range where e ^ 1 we set it to zero and find from 
Eq. ([57)1 the following equation for po^ 



1 



na 
~2 



Po 



(58) 



